Strongly disordered Hubbard model in one dimension: 
spin and orbital infinite randomness and Griffiths phases 
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We study by the strong disorder renormalization group (RG) method the low-energy properties 
of the one-dimensional Hubbard model with random-hopping matrix-elements tmin < t < t,nax, and 
with random on-site Coulomb repulsion terms < Umin < U < Umax. There are two critical phases, 
corresponding to an infinite randomness spin random singlet for strong interactions {Umin > tmax) 
and to an orbital infinite randomness fixed point for vanishing interactions {Umax/tmax — > 0). To 
each critical infinite randomness fixed point is connected a Griffiths phase, the correlation length 
and dynamical exponent of which have well defined asymptotic dependences on the corresponding 
quantum control parameter. The theoretical predictions for the scaling in the vicinity of the critical 
points compare well to numerical RG simulations. 

PACS numbers: 75.10.Nr, 75.50.Lk, 64.60.Ak 
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I. INTRODUCTION 

The Hubbard model is one of the simplest quantum 
many-body models of interacting fcrmions. Detailed ex- 
act information in one-dimension (ID) for the pure model 
is provided by Bethe Ansatai, bosonization^ and by fac- 
torizing the wave- function^. For repulsive interaction 
U > 0, the ID model is the simplest system with a 
Mott transition at half-filling, having a finite charge gap 
Ac > 0. At the same time the spin gap is vanish- 
ing, and the spin-spin correlation function has quasi-long- 
range order. 

Quenched, i.e. time- independent disorder is an un- 
avoidable feature of real materials and there are quasi- 
ID systems, such as tetracyanoquinodimethan (TCNQ) 
compounds^ for which the Hubbard model with ran- 
dom parameters is conjectured to be relevant. Earlier 
theoretical studies^iSii based on the real-space renor- 
malization group (RG) method^, quantum Monte Carlo 
(QMC) simulations^ and density matrix renormalization 
(DMRG)^ are primarily interested in the effect of weak 
random potentials of strength e, in order to understand 
the metal-insulator transition in the presence of interac- 
tions. Repulsive interactions turn out to contrast with 
attractive interactions: like the non interacting model^, 
the system with repulsive interactions is in the insulating 
phase^ whereas a metal-insulator transition is predicted 
for strong enough attractive interaction. For repulsive 
interactions, the charge gap is reduced by the random 
potential and for strong enough e there is a transition 
from a Mott insulator (in which the charge gap is finite 
and the spin-spin correlation length is infinite) to an An- 
derson insulator (with a vanishing charge gap and a finite 
spin-spin correlation length) . 

Other investigationsiSiii of insulating disordered ID 
systems, motivated by the understanding of the NMP- 



TCNQ compound^^, and more recently-'^^ by the inor- 
ganic spin-Peierls compound CuGeOa, have relied on 
the evaluation of the uniform paramagnetic susceptibility 
within strong disorder. The proposed effective model^ is 
a random exchange 5 = 1/2 antiferromagnetic Heisen- 
berg spin chain. The RG introduced by Ma, Dasgupta 
and Hui^iil to study this model was first implemented 
numericalljii^, and then analytically by Fisher—. Ac- 
cording to the RG results, the ground state consists of 
a collection of randomly distributed singlets. Disorder 
grows without limit as the energy is reduced. The result- 
ing singular behavior is typical of an "infinite random- 
ness" fixed point: the low-temperature uniform suscepti- 
bility diverges asi^: 
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which has been first found in numerical RG studiesi^. 
On the other hand the experimentally measured suscep- 
tibility in TCNQ materials is accurately described hyi^ 
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with l3 w 0.55-0.9. The results of QMC simulationa^ are 
in qualitative agreement with the power-law behavior in 
Eq. lI2Jl. A similar behavior was found more recently in 
the insulator o-TaSs where the power-law behavior in the 
susceptibility matches accuratelyiS, both the power-law 
specific heat and the power-law electronic spin resonance 
spectroscopy signalii. 

It is of importance to understand the low-energy prop- 
erties of the Hubbard model in the presence of differ- 
ent types of randomness and arbitrary interactions. As 
shown in Ref. |.18| dilute impurities in an interacting 
quasi- ID conductor can stabilize a "bounded" Luttinger 
liquid on ballistic finite size segments in between two im- 
purities that constitute infinite barriers at low energjiiS. 
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We consider in the following a similar problem in the 
commensurate case and in the regime of strong disorder, 
likely to be relevant to a finite concentration of impu- 
rities. The combined effects of disorder in the hopping 
matrix-element and in the interaction result in spin and 
orbital random singlets, and spin and orbital Griffiths 
phases. The transitions and cross-overs between them is 
the subject of our article. For completeness we also con- 
sider the possibility of diagonal disorder in the form of a 
random potential, giving rise to an Anderson insulator. 

The structure of the article is the following. The ran- 
dom Hubbard model is introduced in Sec. ^] The basic 
steps of the strong disorder RG method together with 
the overall phase diagram of the model are presented in 
Sec. mil Results of numerical renormalization group cal- 
culations are given in Sec. IIVI and discussed in Sec. 



III. STRONG DISORDER RG, AND OVERALL 
PHASE DIAGRAM 

To analyze the low-energy properties of the ID random 
Hubbard model we use the strong disorder RG methodic, 
in which we first select the interaction term (being the on- 
site Coulomb interaction, the hopping integral or the on- 
site potential) with the largest parameter in the Hamil- 
tonian in Eq. that defines the energy scale fl. The 
largest coupling with strength is eliminated, which re- 
sults in new parameters between the remaining degrees 
of freedom, that are calculated in perturbation. We first 
set emax = and thus omit the random potential. 

A. Interaction dominated region 



II. MODEL 

We start with the Hamiltonian of the ID Hubbard 
model: 

L 

*=i <y=hl 

+ ^ U^Ui^^n.^i + ^ CiC+^Q,^ + H' (3) 



We start with a disorder in which the interaction plays 
the dominant role, so that even the smallest interaction 
is larger than the possible maximal value of the hopping 
integral (see Fig. ^) : 

Urain ^ ^max • (5) 

In this case only interactions are decimated out 
in the first stage of the RG, corresponding to U- 
transformations. 



1. U -transformation 



in which the extra term Ti! is generated during renor- 
malization, so that 7Y' = in the initial condition of the 
physical model. The operators c^^ and Ci^o- create and 
annihilate a spin-cr fermion at site i = 1,2, ... ,L, and 
ni^u = cf^Ci^rj. We restrict ourselves to repulsive on-site 
Coulomb interactions being independent random num- 
bers such that < Umm < Ui < Umax, choscu in the 
distribution 



V(U) = Ta{U, Umin, Umax) = " 



{Ur, 



U„ 



(4) 



Similarly, the exchange integrals tmin < ti < tmax 
are independent and identically distributed random 
numbers, and we consider the distribution 'P{t) — 
J^at{t, tmin, tmax)- We havc also included in the Hamil- 
tonian in Eq. ijSJl a symmetrically distributed random 
potential, —Cmax < e < emax- The absolute value |e| 
of the diagonal disorder potential is drawn in the dis- 
tribution V{\e\) = (|e|, 0, Emaa;)- Wc usc the Same 
exponent, a = at = = ajj for the hopping, diagonal 
disorder potential and interaction distributions. The pa- 
rameter = var[ln?7] = var[lnt] = var[ln |e|] measures 
the strength of disorder for C/„iin = tmin = 0, (var[x] 
stands for the variance of a;). A uniform distribution cor- 
responds to a = 1 . We restrict ourselves to the half- filled 



Let us consider the first decimation steps, when the 
largest interaction term is, say, U2 — the hopping in- 
tegral between sites 1 and 2 is denoted by ^i, and it is 
denoted by t2 between sites 2 and 3. The double oc- 
cupancy of site 2 is forbidden (at least up to order of 
0(max(ii, i2)/f/2)) because of the large value of U2- The 
Hamiltonian is then projected on the subspace without 
double occupancy. A virtual exchange interaction is then 
generated between the spin at site 2 and between the 
neighboring sites. If site 1 contains two fermions the ex- 
change interaction is given in a second-order perturbation 
calculation as: 

2tl 
U2 ' 



Ti'i = ^ J1S1.S2, Ji 

Si 



(6) 



L. 



and similarly between sites 2 and 3. New antiferromag- 
netic Heisenberg exchange interactions are thus gener- 
ated during renormalization. All the interactions are first 
decimated out with the condition in Eq. ((SJl. If we per- 
form the U transformation at both ends of a link we 
obtain for the final exchange term: Ji « 2t\/U2 + 2t\/Ui 
and in Ti'i in Eq. jnj there is no sum over Si, since dou- 
ble occupancy of site 1 is also forbidden. Thus after elim- 
inating all the on-site Coulomb interactions we are left 
with the Hamiltonian of the random antiferromagnetic 
Heisenberg chain: TL' — X^i^i Ji^i-^i+i- This model has 
been thoroughly studied by strong disorder RG^"*. The 
corresponding RG-step is called as J-transformation. 
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(Color online) Schematic representation of the different hopping and on-site interaction energy distributions (drawn 
I), the corresponding phases and the transition or cross-over between them. 



2. J -trans formation 

In the decimation procedure the two spins coupled by 
the strongest exchange coupling, say with J2 — fi, form a 
singlet and are eliminated, while between the remaining 
sites, 1 and 4, a new exchange coupling is generated, 
given by J « Ji J3/2J2 , where Ji and J3 are the original 
couplings joining to J2. 

3. Random singlet spin infinite randomness phase 

The coupling distribution broadens without limit as 
the renormalization goes on and the energy scale is low- 
ered, corresponding to a spin infinite randomness fixed 
point for the distributions on Fig.^. At this fixed point, 
the length-scale L and the energy scale A (the smallest 
gap) are related by 

-lnA-Li/2, (7) 

and the distribution function -Pl(A) scales like 

ln[Li/2p^(A)] ^/(L-i/2inA) . (8) 

The ground state is the so called random singlet phase, 
consisting of effective singlets between sites arbitrarily far 
apart. The average spin-spin correlation function decays 
like C(r) ~ r~^. The low-temperature susceptibility has 
a Curie-like form with logarithmic corrections, as given 
in Eq. The predictions of the asymptotically exact 



RG have been confirmed by numerical simulations 
see however Ref. [23|. 



4- Random dimer phase 

To introduce the notion of Griffiths phasc^'' that will 
be useful for discussing the random Hubbard model, we 
consider now an enforced dimerization, with the cou- 
plings Jo at odd positions, and Jg at even positions taken 
from different distributions. The control parameter Sdim 
is defined as ddtm = [lnJo]av - [lnJe]av , where [. . .]av 
stands for an average over quenched disorder. This type 
of model is obtained through renormalization of the ran- 
dom Hubbard model if the hopping terms and/or the 
interactions in Eq. (j3Jl are dimerized and if the relation 
in Eq. Q is satisfied. The dimerized (non-random) Hub- 
bard model is used to describe the low-energy properties 
of Bechgaard salts^^ and other quasi-lD systems^. 

The dimerized random Heisenberg chain is renormal- 
ized to a zero energy fixed pointld*22,, with short range 
spatial correlations i.e. the correlation length is ^ < 00. 
This is the so-called random dimer phase, being a Grif- 
fiths phase consisting of long singlet bonds involving the 
Jo couplings and short singlet bonds involving the Jg cou- 
plings. The average size ^ of the latter is used to define 
the average correlation length, which behaves for small 
dimerization as 

^-5,1. (9) 
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On the other hand, the length of the long singlet bonds 
is proportional to the size L of the system, and the low- 
energy excitations goes to zero as 



A 



(10) 



Here 2 < cx3 is the dynamical exponent which depends 
on the value of Sdim, thus on the distribution of disor- 
der. Close to the random singlet phase, \5dim\ ^ 1, the 
dynamical exponent diverges as 

^ ~ ^dL ■ (11) 

If the low-energy excitations are localized, which is the 
case in the random dimer phase, the gap distribution 
Pi (A) scales like 



Pi(A)=i^5(L^A)^A 



(12) 



The low-temperature uniform susceptibility in the ran- 
dom dimer phase behaves like x(r) ^ i/yi^i/^ ^ which 
is in the same form as the experimental result in Eq. ^ , 
with /? = 1 — 1/z. In the random singlet phase, Eas. (|10|) 
and ifT^ are replaced by Eqs.Q and ^ respectively. 



B. Hopping dominated region 

Some hopping matrix elements are also decimated as 
51 is lowered below tmax, if the condition set by Eq. (O 
on the strength of the on-site Coulomb interaction is not 
satisfied, which we consider now. 



However, 
order: 



if Ji = and/or J3 = we obtain to third 



12 



(16) 



The ^-transformation has an influence on the renormal- 
ized value of the other parameters, in particular on the 
on-site Coulomb interaction. In particular, both t- and 
f7-transformations are generated at some stage of the 
renormalization for any finite value of Umax > 0. 



2. Orbital infinite randomness fixed point: 



Umax — emax — 0. 



The problem reduces to the random tight- 
binding model and the RG process involves solely 
t-transformations if Umax — (and tmax = 0). The fixed 
point of the RG is now an orbital infinite randomness 
fixed pointSSiSS, which is isomorph to the fixed point 
of the random antiferromagnetic Heisenberg chain, as 
discussed in Sec. IIII A 31 In particular the relations in 
Eqs.Q and © remain valid. The ground state of the 
system is made of i- frozen pairs of fermions, where the 
length of a pair can be arbitrarily large. 



1, t -transformation 



C. Interplay between interaction and hopping 



If the energy corresponds to a hopping term, say 
n = t2, then two fermions (one with cr =t and the other 
with a =i) are localized on the bond between sites 2 and 
3, and these sites are eliminated. A small hopping term is 
generated between the remaining sites 1 and 4, which is 
given in second order perturbation by < w —^1^3/^2 • This 
transformation has an effect on the value of the interac- 
tion term, say at site 1, provided Ui has not yet been 
decimated out. For its renormalized value we obtain: 

t2 



(13) 



Similarly the random potential at site 1 is modified and 
transformed as 



ei « ei 



t2 



(14) 



Formula analogous to Eas. ll3|l and p4|l hold also at site 
4. 

Finally, if two J-couplings (denoted by Ji and J3) have 
already been generated at sites 1 and 4, an effective ex- 
change interaction is also generated between sites 1 and 
4: 



J : 



J1J3 
16|i2 



(15) 



In the general situation Umm < tmax, Umax > where 
the t- and [/-distributions overlap (this corresponds to 
Figsn]b,c,d and e), both U and t terms are decimated 
during renormalization and J-transformations are also 
carried out. The ground state of the system is a mixture 
of J-singlets and frozen i-pairs and the average size 
of the singlets, and that of the i-pairs are finite. The 
largest of the two defines the correlation length. At the 
same time the energy scale decreases to zero with the size 
of the system as in Eq. H10|l . We are thus in a Griffiths 
phase, analogous to the random dimer phase described in 
Sec. IIII A 41 which can be divided in two regions. The on- 
site Coulomb interaction plays the dominant role in the 
"spin Griffiths phase" corresponding to £,j > for the 
coupling distributions on Fig. and d, whereas hopping 
is dominant in the "orbital Griffiths phase" , correspond- 
ing to > for the coupling distributions on Fig. ^ 
and e. There is a cross-over, but no sharp transition 
in between these two Griffiths phases, similarly to that 
observed in other random quantum models, such as the 
random dimerized S = 1 chain'^'^ or random Heisenberg 
ladders^. In the following subsections we analyze the 
properties of the Griffiths phases in the vicinity of the 
infinite randomness fixed points. Numerical results far 
from the random critical points are presented afterwards. 
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1. Spin Griffiths phase: < t,nax — Umi,i < Umax 

Let us start to analyze the properties of the system in 
the vicinity of the random singlet fixed point, when < 
imax — U,nin ^ Umax- We define the control-parameter 
du as the fraction of non-decimated /7-terms at = tmax- 

5u = P{U)dU = f ^"""'^""" V . (17) 

Here the second relation holds for the power-law distri- 
bution in Eq. 

As we discussed in Sec. IIII Bl some i-terms are also 
decimated out as fl is lowered below tmax, and the den- 
sity pt of the frozen i-terms is given by pt ~ Su- The 
typical length-scale in the system is given by the typical 
distance between two frozen t-terms, and is thus 

^typ ^ I/Pt ^ S^' . (18) 

On the other hand the average correlation length is given 
by the average distance between spins forming a singlet 
through a J-coupling. In the random singlet phase, i.e. 
for Sjj — 0, the singlets are formed either from odd or 
from even bonds. For S > the correlated segments are 
broken if two neighboring f-terms are frozen, which leads 
to a change of the parity of the singlets, see the discussion 
in Sec. IIII A 41 Thus we obtain for the average correlation 
length 

^ ~ 1/p? ^ 6^' . (19) 

Finally, one considers the typical value of the gap, which 
is related to the value of the gap in the random singlet 
phase of size ~ Eqs. {Tjl, and lead to 

z^S^K (20) 

Now comparing Eas. lfT^ and ^^U\i with those of Eqs.Q 
and Hll|l in the random dimer phase we conclude that Su 
in Eq. Ijl?!) plays the role of the control parameter in the 
spin Griffiths phase. 

2. Orbital Griffiths phase - < Umax ^ tmax 

Second, we consider the behavior of the system in the 
vicinity of the orbital infinite randomness fixed point, 
< Umax ^ tmax- In the initial steps of the renormaliza- 
tion solely t-terms are decimated out, until the energy- 
scale is lowered below Umax-, when also {/-terms are 
eliminated. The density of decimated U sites is given 
by pu ^ I In t/maa; I which foUows from Eq. O, with 
the correspondences A ~ Umax and L ~ p^^. We have 
thus for the correlation length (the average size of the 
t-singlets): £,t ^ p^^ ^ \ \nUmax\'^- Comparing to Eq. 
@ in the random dimer phase we identify the control 
parameter in the orbital Griffiths phase as: 

St = \HUmax/tmax)\~^ - (21) 



Repeating the arguments used in the vicinity of the 
random singlet phase we obtain that the relations in 
Eas. H18|l . H19|l and (|20|1 remain valid by simply replac- 
ing Su by St- 



D. Role of the random potential - Anderson vs. 
Mott insulator phases 

In order to determine the overall phase diagram of the 
random Hubbard model we consider also the role of the 
random potential {ei} which plays a special role in the 
renormalization in ID. If at some stage of the renor- 
malization one random potential term is decimated, say 
£2 = 0, then two fermions are frozen at site 2, and this 
site is eliminated. This transformation does not influ- 
ence the value of the parameters at the neighboring sites 
because no hopping is generated between the remaining 
sites 1 and 3. Transport is thus blocked at this site. The 
random potentials are just slightly modified by other t- 
transformations (see in Eq. H14|l 'l. The random poten- 
tials are thus also decimated at some stage of the RG if 
<^max > Umax > 0. As a consequence a finite fraction of 
sites is frozen and the system scales into the Anderson 
insulator phase. In the other limiting case, Cmax < Umin, 
the behavior of the system is dominated by the on-site 
Coulomb repulsion, thus no fermions are frozen due to 
the random potential and the system scales into the Mott 
insulator phase. 



E. Phase diagram 

We conclude this section by presenting the phase dia- 
gram of the random Hubbard model in ID, which con- 
stitutes the main result of this paper. First we con- 
sider in Fig. the effect of a random potential, when 
the disorder in the other terms is weak, so that {tmax — 

tmin) / {tmax ~t~ ^min) ^ {Umax Umiri) / {Umax ~t- Umin) — 

Dt,u ^ ^max- Previous numerical workSii considered 
non-random U and t, thus Dt^u = 0. Here we extend 
these results to weak disorder. The random potential 
plays a dominant role for tmax /Umin > 1, when the sys- 
tem is in the Anderson insulator phase, with a finite 
spin-spin correlation length ^ < oo and with a vanish- 
ing charge gap Ac = 0. If Dt^u = and Smax/Umin < 1 
the system is in the Mott insulator phase, ^ = oo but 
with Ac > 0. We expect this scenario to remain valid 
even for weak disorder {Dt.u ^ ^max) provided hopping 
dominates over Coulomb repulsion, thus tmax /Umin > 1- 
If, however the on-site Coulomb repulsion is sufficiently 
strong {tmax/Umin < 1 and tmax/Umin < 1 ) the sy stcm 
is in the spin random singlet phase, see Sec. IIII A 31 with 
a divergent spin-spin correlation length, — oo and with 
a vanishing charge gap, Ac — 0. Transition between the 
Mott and the Anderson phases is found to be controlled 
by a conventional random fixed point-''', whereas tran- 
sitions from the spin random singlet phase are likely of 
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FIG. 2: (Color online) Schematic phase diagram of the ran- 
dom Hubbard chain, a) Including a random potential and 
with weak disorder in U and t. b) Without random potential 
and with strong disorder in U and t. This latter corresponds 
to the numerical work in the paper. OIRFP stands for orbital 
infinite randomness, OG for orbital Griffiths (for the distribu- 
tions on Fig. SG for spin Griffiths (for the distributions 
on Fig. SID for spin infinite randomness (for the distri- 
butions on Fig. ^). See the text for the properties of the 
different phases. We suppose as in Fig. Q coupling distribu- 
tions with constant values of Umax ~ Umin and tmax ~ tmin 

while Umax and tmax are varying. 



infinite randomness type. 

In Figl^D we present the phase diagram without a 
random potential, emax = 0, but for sufficiently strong 
disorder in U and i as a function of r = Umax /tmax- 
For Umax /tmax = the systcm is in the orbital infinite 
randomness fixed point, the properties of which are de- 
scribed in Sec. IIII B 21 In the other limiting case with 

Umtn > tmax, which COrrCSpOuds to Umax/tmax > > 0, 

the system is in the spin random singlet phase (see in Sec. 
IIII A 3p . For < Umax/tmax < fc the system is in the 
Griffiths phase, which is divided in an orbital Griffiths 
phase (Sec. IIII G 2|) for < Umax/tmax <C Tc, and in a 
spin Griffiths phase (Sec. IIII G 1|) for <^ Umax/tmax < 
Tc- In the following section we study now numerically the 
properties of the phase diagram in Fig[5l3. 
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FIG. 3: (Color online) The fractions of different RG transfor- 
mations as a function of the log-energy scale, P = — Infi, for 
the random ID chain. 



tained similar results for smaller values of a. The length 
of the chain is varied up to L = 2048 and the decima- 
tion is performed up to the last remaining particle or 
spin singlet. The gap A at the last step of renormaliza- 
tion is identified as the gap of the random chain and we 
have also evaluated the length scales and O from the 
density of the decimated t and J terms respectively. 



A. Qualitative features of the RG flow 

To illustrate the behavior of the RG flow we have cal- 
culated the fraction nt{T)5T of ^-transformations, the 
fraction nii{^)5T of [/-transformations, and the fraction 
nj{T)5T of J-transformations in an interval [r,r -I- JF], 
where F = — In f2 is the log-energy scale. The quantities 
are normalized in such a way as ntiV) -|-n[/(F) ~\-nj{T) = 
1. 

The variations of nt{T), nu{T)^ "-./(r) for the ID chain 
are shown on Figs.|21for Umax =0.1, tmax — 1 and a = 1 



with {/„ 



tr^ 



0. As seen in this figure only t- 



terms are decimated out at the beginning of the RG, 
and the decimation of the U -terms starts only when Vt is 
lowered below Umax- J-terms are also generated through 
[/-decimations, which are then also decimated. The U- 
and i-terms gradually die out by further decreasing 
and one is left with the J-transformations at low energy. 



B. Properties of the Griffiths phases 



IV. NUMERICAL RESULTS 



1. Dynamical exponent 



We consider typically 250000 independent random 
samples taken from the distribution in Eq. to evaluate 
the dynamical exponent in the numerical implementation 
of the strong disorder RG method. We use a smaller 
statistics of ^ 10000 samples for the evaluation of the 
correlation lengths. The numerical results are presented 
only for the uniform distribution with a = 1, but we ob- 



The dynamical exponent z is determined from the dis- 
tribution of the gaps, which in the Griffiths phases follows 
the scaling form given in Eq. (|12|l . For a finite system 
this procedure results in effective L-dependent dynamical 
exponents, that are then extrapolated to L = cxd. The 
calculated effective dynamical exponents for the power- 
low distribution in Eq. Q with tmin = Umin = and 
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FIG. 4: (Color online) Dynamical exponent as a function of 

r — Umax/tmax for different system sizes {Umin = tmin — 0, 
a = 1). 



FIG. 5; (Color online) Scaling plot of the dynamical exponent 
in FigHlin the spin Griffiths phase, with 5u in Eq. If 711 . 



a = 1 for various values of r = Umax/tmax are shown 
in Fig. 21 For these distributions the infinite random- 
ness fixed points are located at r = (orbital infinite 
randomness fixed point), and at r = oo (spin random 
singlet, infinite randomness phase). The dynamical ex- 
ponent z is formally infinite at both fixed points, which is 
compatible with the fact that the effective z is increasing 
with L with no sign of saturation. It is also interesting 
to note that the maxima of the effective exponent in the 
orbital Griffiths phase for a size L are approximately the 
same as the maxima of the same curve in the spin Grif- 
fiths phase, however with a size L/2. This is because the 
spin Griffiths phase involves approximately the double of 
RG steps (in the spin Griffiths phase one makes also a 
set of extra {/-transformations). 

To have a qualitative analysis of the data we recall that 
close to the infinite randomness fixed points the effective 
exponent z{5, L) obeys the scaling form: 

z{5,L)^r^z{L^'^5) , (22) 

in which we have incorporated relations in Eq. 1111) and 
in Eq. Furthermore, the scaling function behaves 

like z{x) ~ X for small x, leading to an ^-independent ~ 
L^/'^ effective dynamical exponent at the critical points. 
The scaling form in Eq. (|22|) is well verified in the vicinity 
of the two infinite randomness fixed points, keeping in 
mind that the control parameter 5 is given in Eas. (|17|) 
and 1)21(1 for the spin and the orbital Griffiths phases 
respectively. We obtain a satisfactory agreement between 
the numerical results and the analytical calculations in 
both regions (see Figs. [S] and (HI. 

2. Correlation length 

The correlation length is finite in the Griffiths phases 
and given either by the average length of the spin 
singlets in the spin Griffiths phase, or by that of the 
orbital singlets in the orbital Griffiths phase. In terms 
of the appropriate control parameter, given either in Eq. 
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FIG. 6: (Color online) The same as in Figl^ in the orbital 
Griffiths phase, with &t in Eq. 12111 . 

(|17|l or in Eq. (|21|) . the correlation lengths are divergent 
in the vicinity of the infinite randomness fixed points, 
as given in Eq. 0. This relation has been verified by 
numerical simulations (see Fig|7|), in which - in order to 
get rid of finite size effects - we used large finite systems 
with length up to i = 2048. For the spin Griffiths-phase 
we used a distribution of the form in Fig^D, whereas 
for the orbital Griffiths phase the distribution is given 
in Fig^. As seen in the figure satisfactory agreement 
with the theoretical prediction is found, at least for small 
enough 5 and for large L. 



V. DISCUSSION AND CONCLUSION 

In this paper we have studied the low-energy, low- 
temperature properties of the one-dimensional Hubbard 
model in the presence of quenched disorder. Contrary 
to previous studies^*^ in which disorder was compara- 
tively weak and realized in the form of a random po- 
tential, we considered here disorder both in the hop- 
ping integrals and in the on-site Coulomb repulsion. Due 
to the sufficiently strong disorder, the system is always 
gapless, which means that the low-energy properties are 
controlled by zero-energy fixed points, for any value of 



8 




(1/5)' 

FIG. 7: (Color online) Correlation length in the spin Criffiths 
phase and the orbital Criffiths phase given by the average 
length of the spin and orbital singlets respectively as a func- 
tion the inverse of the square of the control parameter (see 
Eqs. 1)17^ and H21|l respectively). The theoretical predictions 
correspond to the solid line ^ = . 

the quantum control parameter. In the framework of 
the strong disorder renormalization group that we have 
adapted here to the random Hubbard model, we have 
identified two distinct critical phases, which are con- 
trolled by infinite randomness fixed points. For dominant 
on-site repulsion the system is a spin random singlet in- 
finite randomness phase, whereas for vanishing on-site 
Coulomb interaction the system corresponds to a ran- 
dom tight-binding model the properties of which are con- 
trolled by an orbital infinite randomness fixed point. 

The two infinite randomness critical phases are sep- 
arated by two different Griffiths phases, between which 
there is a smooth cross-over, but no sharp transition. The 
strong disorder RG method is expected to provide asymp- 



totically exact results close to the critical phases, where 
the correlation length is divergent. We have made analyt- 
ical predictions, both for the divergence of the correlation 
length and that of the dynamical exponent. These results 
were compared to large scale numerical RG calculations 
and satisfactory agreement is found. 

To conclude, we mention possible extensions of our 
work. It is of interest to study the weak-to-strong disor- 
der effects in the random Hubbard chain. For example 
the charge gap in the Mott insulator phase is expected to 
be robust against weak disorder, but strong disorder will 
destroy it, as shown by our RG results. It would be also of 
interest to study the combined effects of disorder realized 
at the same time in different parameters (potential, hop- 
ping and interaction). Finally, one may consider exten- 
sions of the model to higher dimensions. A new physics is 
expected already in two dimensions because the interac- 
tion dominated phase that maps to a random Heisenberg 
model has a conventional random fixed point"^^, whereas 
for vanishing on-site Coulomb repulsion the system is in 
a logarithmically infinite randomness fixed point^i. How- 
ever, for a quasi-one dimensional system, made of weakly 
coupled chains, the transverse exchange couplings of the 
low energy Heisenberg model are much smaller than the 
longitudinal ones because of Eq. The one-dimensional 
spin Griffiths phase due to strong on-site Coulomb inter- 
actions is thus expected also in the quasi-one-dimensional 
system with sufficiently small interchain couplings. 
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